Strong coupling of a qubit to shot noise 
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We perform a nonperturbative analysis of a charge qubit in a double quantum dot structure 
coupled to its detector. We show that strong detector-dot interaction tends to slow down and halt 
coherent oscillations. The transitions to a classical and a low-temperature quantum overdamping 
(Zeno) regime are studied. In the latter, the physics of the dissipative phase transition competes 
with the effective shot noise. 
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The study of fluctuations and noise provide deep in- 
sights into quantum processes in systems with many de- 
grees of freedom. If coupled to a few-level system such 
as a qubit, fluctuations usually lead to destabilization of 
general qubit states and induce decoherence and energy 
relaxation. One important manifestation is the back- 
action of detection on qubits 0. This topic has been 
extensively studied in the regime of weak coupling be- 
tween qubit and noise source |2j ■ It has been shown that 
the qubit dephases into a mixture of qubit eigenstates 
(dephasing), whose classical probabilities thermalize to 
the noise temperature at a longer time scale. Meso- 
scopic charge detectors such as quantum point contacts 
(QPCs) 0] and radio-frequency single electron transistors 
(rf-SETs) S, whose low-temperature noise is shot noise 
0, E| , are particular powerful detectors as they provide 
high resolution and potentially reach the quantum 
limit. A particular attractive regime for qubit applica- 
tions is the QND regime, realized if the qubit Hamilto- 
nian and the qubit-detector coupling commute 0, . 

We study a quantum point contact potentially strongly 
coupled to the coordinate (left / right) of a double quan- 
tum dot charge qubit 0, by a nonperturbative ap- 
proach involving the Gaussian and nonintcracting blip 
approximations. We analyze the qubit at the charge de- 
generacy point, where the two lowest energy eigenstates 
are delocalized between the qubits. In the weak coupling 
regime, low-temperature relaxation would thus always 
delocalize charge. We show that, in strong coupling, the 
qubit state gets localized in one of the dots. Localization 
is manifest by a suppression of both the coherent oscilla- 
tions and the incoherent tunneling rate. This "freezing" 
of the state also applies a high bias and can e.g. lock an 
excited state. Thus, in the strong coupling regime, the 
environment naturally pushes the physics to the QND 
limit even if the bare Hamiltonian is not QND. We point 
out the analogy of this physics to the case of the dissipa- 
tive phase transition in oscillator bath models , which 
in the QPC competes with the nonequilibrium induced 
by the voltage driving the shot noise. 

We consider the case of a degenerate two-state system 
(TSS), realized by the charge states in a double quantum 
dot structure (see Figure Q). These charge can be read 
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FIG. 1: Schematic view of the double dot system analyzed 
see e.g. Refs. 0,0- The QPC and rf-SET detectors can be 
used alternatively, both options are discussed in the paper. 



out by the current through a nearby quantum point con- 
tact. The Hamiltonian for the TSS with time-dependent 
fluctuation e(t) reads 
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(1) 

In the last step of eq. JTJ, we applied a Polaron trans- 
formation |12| introducing the fluctuating phase <p = 
/ dt'e(t'), with i(t) — e + Se(t), for the tunneling matrix 
elements in the qubit. The microscopic foundation of the 
noise term 5e(t) for a QPC is give n in Refs. an d 
for an SET in Refs. [Hill mi!! 

Without loss of generality, we assume (ct z (0)) = 1. We 
can now formally solve the Liouville equation. The ex- 
pectation value of <7 2 , the difference of occupation prob- 
abilities of the dots, satisfies a closed equation 
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where the second line of eq. (J2J has been derived by as- 
suming that the noise represented by J(t — t') is station- 
ary. This procedure is analogous to the nonintcracting 
blip approximation (NIBA) of the path-integral solution 
of the Spin-Boson model 0, 0] . This automatically in- 
cludes a Gaussian approximation to the shot noise 0. 
This approach is nonperturbative in (f> and a good ap- 
proximation in the two cases e = and |ei 3> |A|. 
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We start with the charge-degeneracy case e — 0. Here, 
we can solve eq. (|2l in Laplace space and find 
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with the Laplace-transformed self-energy S(s) = 

oc 

A 2 J dte~ st e Jtyt \ The phase correlation function J(t) as 

o _ 
seen by the dots reads |5j 

oc 

27T /■ , |Z( W )| 2 



J(t) 



KR 



K 



du>- 



(4) 



where Si(oj) is the full current noise in the QPC that for 
sufficient environmental impedance is given [5j by 
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and the transimpedance Z(u>) between qubit and point 
contact. In eq. (J5J, V is the bias voltage of the QPC, Rk 
is the quantum resistance, and D m is the transmission 
eigenvalue of the mth conductance channel. 

Semiclassical limit: We now discuss the resulting dy- 
namics in a number of limiting cases. We start by 
first taking the limit lu — > 0. This corresponds to 
TiA,he -C eVjksT, i.e. the qubit probes the shot 
noise at energy scales much lower than its internal ones. 
Here, the noise expression [cq. (JjJJ] becomes frequency 
independent 0. We can then compute the semiclas- 
sical spectral function J c (t) — —j c t- Here, we have 
assumed a frequency-independent transimpedance con- 
trolled by a dimensionless parameter k, \Z(u>)\ 2 « k 2 R k 

N 

and 7c = 2it 2 k 2 R k S i {0) with 5/(0) = ^E-°m(l - 



D m )eVcoth^y 
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T~ S § ■ The self-energy is then 
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readily calculated and analytical, so we can go back from 
Laplace to real time and obtain 
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We observe that the coher- 



where w cff , c = y A 2 
ent oscillations of the qubit decay on a scale 7~ 1 and 
get slowed down. At j c = 2 A, the damping becomes 
critical and the oscillations disappear, ending up with a 
purely exponential overdamped regime at j c > 2A. This 
crossover corresponds to the classical overdamping of a 
harmonic oscillator. Even in the overdamped regime, the 



qubit decays exponentially to (<r z (t)) — ► at long times, 
e.g. it gets completely mixed by the shot noise, whose 
noise temperature is high kBT nomc ~ max{ey, kgT} ^> 
fiA. Note that it is possible to discuss the overdamped 
regime, where 7 C is not a small parameter and our theory 
is also non-Markovian, see eq. J2J), capturing the neces- 
sary time-correlations arising in strong coupling. 
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FIG. 2: Semiclassical limit: expectation value (a z (t)) as a 
function of time and varied bias voltages V at e = 0. The 
other parameters are T = 0.1 K, D = 0.1, A = 1.524- 10 9 1/s, 
and k — 0.02. Inset: as a function of time and the QPC 
transmission D with fixed QPC bias voltage eV = 100 hA. 

Figure [21 shows the resulting dynamics in the one- 
channel case. With increasing bias voltage V over the 
QPC, the expectation value (& z (t)} drops down quite 
fast. The transmission D of the QPC has also an impor- 
tant impact on the stability of the oscillations of (a z (t)) 
(see inset of Figure EJ|. At D — 0.5, the expression for 
5/(0) has a maximum, therefore the oscillations are there 
maximally suppressed. 5/(0) represents the shot noise of 
the QPC in the low frequency regime ^^|. The more 
noise the QPC provides, the quicker the oscillations de- 
cay. Note that changes in the QPC transmissions (and 
therefore the Fano factor) do not play any role other than 
entering the total noise level. 

Quantum limit: Now, we let T — > and leave u arbi- 
trary. 5/(w) reads in this limit 
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This shape is dominated by two terms, which resem- 
ble the Ohmic spectrum at low T, Sq cx oj0{uj) with 
shifted origins of energy. For computing the quantum 
correlation function J q (t), an ultraviolet cutoff lj c has 
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to be introduced, which physically originates either from 
the finite bandwidth of the electronic bands in the mi- 
croscopic Hamiltonian or from the high-frequency limi- 
tations of the transimpedance Z(u>). We end up with the 
long-time limit for J q (t) applicable at hA <C eV 



Tl I UJ r , 



-j q t + ia 3 . (8) 



This holds for any number of channels, for simplicity we 
concentrate on the single-channel case with a Fano factor 
then is given by F = 1 — D, which we use from now on. 
Here, we can introduce o?2 = g = 16itk 2 D, the dimen- 
sionless conductance as seen by the qubit, a\ = g"fD, 
a 3 = 7r<7 / 2 and 7 9 = irg(l — D)eV/2K. The resulting 
self-energy is now non-analytical 
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where we have introduced the effective tunnel splitting 
A 2 ff = A 2 e~" lgD (^2) ?(-gD+l). In our regime, 
LD C 3> eV/h ^> l/t ~ A, this expression resembles the 
renormalized A of the Spin-Boson model |l8| and we 
have A e ff <C A. This is a sign of massive entanglement 
between system and detector j^, HJ. Note that simi- 
lar to the adiabatic scaling treatment in Ref. 0], the 
NIBA is compatible with forming entangled states be- 
tween system and bath. This has been numerically con- 
firmed, for the Spin-Boson model, in Ref. [2(J. An elegant 
approach to this system reflecting entanglement and use 
of the measurement result in the perturbative regime has 
been given in Ref. j^. The main difference in our shot 
noise case is that the infrared cutoff entering the renor- 
malization and controlling the final expressions appears 
to be V instead of A. In particular, A e g grows with eV, 
which indicates that nonequilibrium shot noise competes 
with the Spin-Boson-like suppression. 

The self-energy is analytical only at F = 1, which 
corresponds to the no-noise case D = 0. Due to the 
generally non-analytic self-energy, it is difficult to com- 
pute the full real-time dynamics by back-transformation 
to the time domain. The structure of the result will be 
(a z (t)) = P cut {t) + P coh (t) + P incoh (t) [3. For our case 
of e = 0, there is no incoherent exponential decay -Pmcoh- 
P cu t is a nonexponential branch cut contribution. In the 
following, we concentrate on the coherent part P co h{t), 
given through the poles Si — — j c g ztiucs of 5 with finite 
imaginary part, and hence this leads to damped harmonic 
oscillations with frequency w e ff and decay rate j c s ■ 

Close to D = 0, we can characterize these poles per- 
turbatively. We find a renormalized oscillation frequency 



Wrff, 



namely cu c ff — Re I y A 2 (l + ^g) — f- J whereas 
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Here, A 2 is de- 



fined as A 2 = A 2 (l + gln^j^jy For arbitrary F 

or _D, we can solve the pole equation numerically, see 
Fig. [2 With the numerical results from Figure^! one can 
again calculate the Laplace back-transformation, where 



the two residues of the kind a_i 



, f (^+7,) 
S ,(2- 9 D)+ 7 , 



have to 



be summed up. This leads finally again to decaying os- 
cillations as already mentioned above. 

We see that at sufficiently strong coupling to the de- 
tector, a finite Fano factor can lead to a complete sup- 
pression of the coherent oscillations, whereas the decay 
rate increases. Both these tendencies together show that 
a finite Fano factor brings the system closer to charge 
localization. In fact, for sufficient damping, we can tune 




0.4 0.6 
QPC transmission D 

FIG. 3: Quantum limit: imaginary parts of the numerically 
determined poles as a function of the QPC transmission D. 
The other parameters are eV = 10 2 hA, Lo c = 10 12 A. Inset: 
real parts of the poles as a function of the QPC transmission 
D. 



the tunneling frequency all the way to zero by increasing 
D. On the other hand, also 7 e g can become very small — 
in these points the detector completely localizes the par- 
ticle up to nonexponential contributions. At other values 
of D, unlike the dissipative phase transition in the Spin- 
Boson model, the hot electrons driving the shot noise 
again drive the relaxation rate close to its bare value, 
and thus this resembles the classical overdamping case. 

This scenario is not limited to e = 0. NIBA permits 
to reliably study the op pos ite regime £ > A as well. As 
already shown in Refs. [ill E^ |. the resulting dynamics 
is dominated by incoherent exponential relaxation dom- 
inating over P co h and P cu t- The relaxation rate is 



T r = 2Re [S(«e + 0)] 



2A 2 ff Re 



(is + 7,) 



gD-l 



(eV) 



gD 
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(10) 

This again demonstrates the slowdown (through A c g ) of 
the decay to the other dot due to the interaction with the 
detector. Notably, this rate does not display standard 
detailed balance at T — 0, rather, around e = 0, the rate 
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is smeared out on a scale of j q , reflecting the role of the 
nonequilibrium shot noise temperature. We have plotted 
this result in Figure 0] 
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FIG. 4: Quantum limit: relaxation rate T r as a function of the 
QPC transmission D. The other parameters are e — 10 A, 
A = 1.524 ■ 10 9 1/s, eV = 10 2 TiA, lj c = 10 12 A. Inset: 
relaxation rate F r as a function of the qubit bias e. Other 
parameters as above, but with D = 0.5. 

Another view on this is that the effective size of the 
noncommuting term between qubit and detector, given 
by A c ff, is reduced, hence the strong interaction brings 
the effective Hamiltonian closer to a QND situation. 

On the other hand, such dynamics is known as 
the quantum Zeno effect. Note that unlike standard 
derivationsdHi^l, this has been derived in a nonpertur- 
bative way, which is consistent with the necessary strong 
coupling and which retains the non-Markovian structure. 

Summarizing the QPC results, we can observe that, on 
the one hand, the system shows traces of the physics of 
environment-induced localization, which competes with 
classical overdamping by effectively " hot" electrons at fi- 
nite voltage and somewhat reinforced at finite Fano fac- 
tor. This can be understood as follows: the dissipative 
phase transition occurs when the environmental noise is 
highly asymmetric in frequency and when the full band- 
width plays a role. At high voltage, the asymmetry of 
the shot noise spectrum is reduced In fact, the j q t 
contribution in the correlation function J q (t) resembles 
the finite temperature term in the correlation function 
of the Ohmic Spin-Boson model — both terms originate 
from the zero-frequency part of the noise. 

A similar analysis on back-action by strong coupling 
of a QPC to a quantum device — there an Aharonov- 
Bohm experiment |2^| — has been done in Ref . . That 
work concentrates on a stationary situation and weak 
hopping into the dot, whereas in our case the dots are not 
connected to leads. The inter-dot interaction however is 
strong and we concentrate on the real-time dynamics. 

These results can be extended to shot noise sources 



other than QPCs. In fact, it may today be quite challeng- 
ing to reach K-values high enough, such that slowdown 
and localization can be observed, when the noise source 
has only a few open channels. An attractive alternative 
is given by readout using metallic SETs fabricated on an- 
other sample layer see Fig. ^ In these devices, there 
is a number of rather opaque conductance channels. 

In that case, we use the expression 0,0,0,0] of the 
voltage noise of the SET (only valid for small frequencies) 



S v {w,wi) = 4 



Tp2 

^SET 
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16wf 
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where £?set = 2 c ET * s ^ ne cnar g m g energy of the SET 
and loi — I/e is the tunneling rate through the SET. 
Then the final result for (a z {t)) is again the same as in 
eq.(JHJ). The difference, of course, is that j c is now defined 

2tt 2 k 2 E 2 

as 7 C = hRKe'hu? ' ^ ne ^ UU ina-ntum mechanical analysis 
in the low-temperature regime works along the same lines 
as the QPC case but goes beyond the scope of this Letter. 

We performed a nonperturbative analysis of the quan- 
tum dynamics of a double quantum dot coupled to shot 
noise. We analyze the crossover from under- to over- 
damped oscillations in the classical case. In the quan- 
tum case, we show that at strong coupling the oscillations 
show the same behavior, competing with a critical slow- 
down similar to the dissipative phase transition. This 
can be interpreted as the onset of a Zeno effect. 
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